Final Analysis and Results of the Phase II SIMPLE Dark Matter Search 
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We report the final results of the Phase II SIMPLE measurements, comprising two run stages of 
15 superheated droplet detectors each, the second stage including an improved neutron shielding. 
The analyses includes a refined signal analysis, and revised nucleation efficiency based on reanalysis 
of previously-reported monochromatic neutron irradiations. The combined results yield a contour 
minimum of a v = 5.7 x 10~ 3 pb at 35 GeV/c 2 in the spin-dependent sector of WIMP-proton 
interactions, the most restrictive to date for Mw < 60 GeV/c 2 from a direct search experiment and 
overlapping for the first time results previously obtained only indirectly. In the spin-independent 
sector, a minimum of 4.7 x 10 -6 pb at 35 GeV/c 2 is achieved, with the exclusion contour challenging 
a significant part of the light mass WIMP region of current interest. 



The search for weakly interacting massive particle 
(WIMP) dark matter remains at the forefront of modern 
physics activity. Estimated to comprise ~ 23% of the 
Universe mass, it is the role of direct detection efforts to 
elaborate its nature, and whether its interaction with nu- 
cleons is spin-independent (SI) or spin-dependent (SD). 
SIMPLE (Superheated Instrument for Massive ParticLe 
Experiments) |1[ is a direct search activity using super- 
heated liquid detectors, and one of only a few in the 
international panorama with sensitivity to the WIMP- 
proton sector of the SD phase space. It is operated at 
the 1500 mwe level of the Low Noise Underground Lab- 
oratory (LSBB) in southern France. 

In [l[ , we reported the first results of a two stage Phase 
II measurement, comprising a 14.1 kgd Stage 1 exposure 
of 15 superheated droplet detectors (SDDs) 043| with a 
total active mass of 0.208 kg. We here provide the re- 
sults of the full Phase II measurement, including a 13.67 
kgd Stage 2 exposure of a second 15 SDD set, together 
with improved neutron shielding and a refined analysis 
of the individual detector run signals, sensitivities, and 
nucleation efficiency. 

A SDD consists of a dispersion of superheated liquid 
droplets homogeneously distributed within a gel matrix, 
which may undergo a transition to the gas phase upon 
energy deposition by incident radiation. Two conditions 
are required for the nucleation of the gas phase of the 
superheated droplets @: (i) the energy deposited must 
be greater than a thermodynamic minimum, and (ii) this 
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energy must be deposited within a thermodynamically- 
dcfincd minimum distance (Ar c ) inside the droplet, where 
A is the nucleation parameter and r c = the thermody- 
namic critical bubble radius. Adjustment of the two con- 
ditions results in the necessity of depositions of order > 
150 keV//im for a bubble nucleation, rendering the SDD 
effectively insensitive to the majority of traditional de- 
tector backgrounds (including electrons, 7's and cosmic 
muons) which complicate more conventional dark mat- 
ter search detectors, leaving only a- and neutron-induced 
events. 

The 15 Stage 2 SDDs were fabricated as described in 
[l|, each containing between 11-19 g of C2CIF5 for a total 
active mass of 0.215 kg; an additional, freon-less but oth- 
erwise identical, SDD again served as an acoustic veto. 
These were initially pressurized to 2.00±0.05 bar, and 
installed at the rate of one per day in a 700 liter water 
pool maintained at a bath temperature of 9.0±0.1 °C, 
this time with the DAQ initiated only after the instal- 
lation of each 8 detector set. The instrumentation was 
identical to that of Stage 1; in contrast to Stage 1 how- 
ever, the SDD pressures were allowed to rise with time in 
order to obtain additional information on the measure- 
ment sensitivity. 

Also in contrast to Stage 1, the water pool rested on 
an additional 10 cm of wood and paraffin, and 10 cm 
of polyethylene, with a rebuilt 50-75 cm thick surround- 
ing water shield. As a result of the seasonal increase 
in water circulation within the mountain, the ambient 
radon level increased to ~ 1000 Bq/m 3 ; continued purg- 
ing of the cavern air reduced this to ~ 100 Bq/m 3 , and 
circulation of the pool water in combination with ra- 
dioassays of the detector construction materials, yielded 
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a Stage 2 a-background estimate, including both pro- 
genitor and daughter decays, of 5.72 ± 0.12 (stat) ± 0.29 
(syst) evt/kgd. 

Extensive Monte-Carlo estimates of the expected neu- 
tron background, which accounted for spontaneous fis- 
sion plus decay-induced (a,n) reactions, and included 
the increased below-pool shielding and new materials ra- 
dioassays, yielded a reduced rate of 0.333 ± 0.001 (stat) 
± 0.038 (syst) evt/kgd; recalculation of the Stage 1 
disposition with the new radioassays yielded a revised 
background rate of 0.976 ± 0.004 (stat) ± 0.042 (syst) 
evt/kgd, with the primary contribution being the con- 
crete. For the improved shielding of Stage 2, background 
neutrons originate mainly from the glass detector con- 
tainment and shield water. 

Stage 2 data was obtained between 12 April - 22 July 
2010. The total exposure was 13.67 kgd, from the de- 
tector installation protocol and mechanical failure of 4 
SDDs during the run as a result of overpressuring; no 
weather-induced data losses occurred. 

Analysis of the Stage 2 signals, as per Stage 1, included 
a filtering of the initial data set (1997 events) via a pulse 
validation routine, a cross-correlation of the remaining 
set in time between all SDDs, and coincidence rejection 
as due to local noise events and that a WIMP interacts 
with no more than one of the in-bath detectors [1]; the 
analysis was improved via a new bandpass filter for noise 
suppression. The signal waveform, decay time constant 
and spectral density structure of the remaining 826 sin- 
gle events were next inspected individually. A particle- 
induced nucleation event possesses a characteristic fre- 
quency response, with a time span of a few milliseconds, 
a decay constant of 5-40 ms, and a primary harmonic 
between 0.45-0.75 kHz; these parameters differ signifi- 
cantly from those of gel-associated acoustic backgrounds 
such as trapped N 2 gas, fractures and local acoustic back- 
grounds such as water bubbles The event-by-event 
analysis permits isolation of the particle-induced nucle- 
ation events with an efficiency of better than 97% at 95% 
C.L. 

Figure 1 displays the signal amplitude and frequency 
for each of the identified 41 particle-induced signal events 
in Stage 2. Following calibrations as described in |l|, a 
nuclear recoil discrimination cut for A < 100 mV was 
again imposed with an acceptance of > 97%, yielding a 
total of 2 events for the entire exposure. 

The two bubble nucleation criteria are thermodynamic 
, so that variation of either temperature or pressure 
modifies the recoil threshold energy and thus the SDD 
sensitivity, as seen in Fig. 2 where the expected varia- 
tion in threshold recoil energies (E*thr) of both neutron- 
induced recoils and a's for several operating pressures is 
shown. The a threshold curve shifts to higher tempera- 
tures with increasing pressure. Since the curves depend 
on A Q, comparison of experiment and theoretical pre- 
dictions with varying A confirmed our measurements Q 
of A = 1.40 ±0.05, yielding no a-sensitivity whatsoever 
above 2.30 bar as observed experimentally. This was then 



used in calculating the ion recoil energy curves shown in 
Fig. 2. For pressures < 2.20 ± 0.05 bar, the threshold 
recoil energy at 9°C remains below 9.0 ± 0.3 keV. 
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FIG. 1: scatter plot of the amplitudes and frequency of the 
primary harmonic of each true nucleation event observed over 
the Stage 2 exposure, with the boxed event for pressures < 
2.2 bar. 
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FIG. 2: variation of the recoil and a energy thresholds with 
temperature for the three C2CIF5 constituents at 2.00 (solid), 
2.50 (dash) and 3.00 (dash-dot) bar, with A=1.40. The freon 
constituents are identified for 2.5 bar. The vertical line indi- 
cates the 9°C measurement temperature; the horizontal line, 
a threshold recoil energy of 8 keV. 

As also seen in Fig. 2, the ion recoil threshold curves 
similarly evolve to higher temperatures and energies with 
pressure increase, rendering the SDDs increasingly less 
responsive to the on-detector neutron spectrum: at 2.50 
bar, this reduction is 30%, consistent with the observed 
absence of any low amplitude events above 2.3 bar when 
weighted by the exposure. 

The pressure records of all SDDs were next inspected 
for evolution during the measurement, and correlated 
with the signal records. Data obtained at pressures 
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greater than 2.20 bar were excluded, reducing the Stage 
2 exposure to 6.71 kgd; correlation with the signal 
record yielded 1 recoil event consistent with the esti- 
mated 2.2±0.3 background neutrons. The Stage 1 events 
were similarly pressure-correlated, reducing the expo- 
sure to 13.47 kgd; reanalysis of the recoil signals via a 
Hilbert transform-based demodulation identified 4 with 
exponential decay characteristic of nonuniform impulses 
observed in acoustic background studies associated with 
SDDs in vibrational contact with their support and air 
bubbles from water inflow, reducing the recoil events to 
10, slightly below the estimated 13±0.6 background neu- 
trons. 

The first Stage 1 results resulted in part from a the- 
oretical bubble nucleation efficiency given by ry(E)=l- 
Ethr/Edep |8|. This r\ however represents only a first 
approximation to the statistical nature of the energy de- 
position and its conversion into heat [9]: a detailed re- 
analysis of previous monochromatic (54 and 149 keV) 
neutron irradiation data [l(| , at 1 and 2 bar as a function 
of temperature, yielded a refined efficiency of 77'=l-exp[- 
r(E/E t /j r -l))] with T=4.2 ± 0.3, independent of pressure. 

We show in Fig. 3 the impacts of the Stage 2 and 
reanalyzed Stage 1 results on SD WIMP-proton scatter- 
ing, together with the competitive results of other direct 
IlT{l3| and indirect [H], [IH experiments. The contours 
are calculated using the previous [l[ Feldman- Cousins ap- 
proach [l6j | based on observing n events against a back- 
ground one systematic uncertainty below the estimated 
neutron-generated recoil background, rf with T=3.6, the 
standard isothermal halo and a WIMP scattering rate 
[l8| with zero momentum transfer, spin-dependent cross 
section for elastic scattering. The form factors of 
[l8| have been used for all odd-A nuclei, with the spin 
values of [1^ used for 19 F; for 35 C1 and 37 C1, the spin 
values are from [2(|, while for 13 C they were estimated 
using the odd group approximation. The Stage 2 result 
is seen to nearly equal the revised Stage 1 result with its 
revised minimum of a p = 9.2 x 10~ 3 pb at 35 GeV/c 2 , 
despite half the exposure. 

The above representation neglects the non-negligible 
spin contribution of the neutron sector in 19 F, which 
is captured in a mo del- independent SD formulation (2p| 
with <jsd ~ [ap <S P >+a„ <S„ >] 2 , where a Pj „ are the 
WIMP-proton, neutron coupling strengths, and <S Pi „ > 
are the expectation values of the proton (neutron) group 
spins. In this representation, experiments define a band 
(single nuclei targets) or an ellipse (multi- nuclei target), 
with the allowed area defined by the intersect of the most 
sensitive results At M^y = 50 GeV/c 2 , com- 

bined with neutron-sensitive XENON 10 [2l|, the allowed 
area reduction is better than 2 /3 compared with Ref. [l[ ; 
masses above or below this choice yield slightly increased 
limits for most all experiments. More relevant would 
however be the model-independent results for M\y ~ 10 
GeV/c 2 , unavailable for the majority of experiments. 

The impact of the results in the SI sector is shown 
in Fig. 4 in comparison with results from other leading 
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FIG. 3: various spin-dependent WIMP-proton exclusion con- 
tours for Phase II, together with the leading direct [riT-[l3T] 
and indirect SuperK IceCube [l5| ] search results; shown 
are the Stage 2 result, the reanalyzed Stage 1 result, and a 
merging of the two. The region outlined in grey is favored by 
cMSSM [H. 



search efforts [12j, |13j, 1 211 - 13 If , again calculated with the 
standard isothermal halo and WIMP elastic scattering 
rate of Ref. [l8j using Feldman- Cousins, a Helm nuclear 
form factor, and r/\ Again, the Stage 2 contour is nearly 
equal to the revised Stage 1 contour with its contour min- 
imum of 7.6xl0 -6 pb at 35 GeV/c 2 . Owing to the low 
recoil energy threshold, both results enter the possible 
light mass WIMP region recently suggested by CoGeNT 

[ail and cresst-ii m. 

A straightforward combination of the two results us- 
ing the Feldman-Cousins approach, based on 11 candi- 
dates with an assumed background 1 cr(syst) below the 
expected total background, yields the "merged" contours 
indicated in each of Figs. 3 and 4; in the SI case, the 
contour minimum drops to 4.7xl0~ 6 pb and the result 
is in tension with the recent reports of CoGeNT [3lj |. 
DAM A/LIBRA [H and CRESST || regarding light 
mass WIMPS, using a significantly different technique 
with different systematics than the XENON [32| and 
CDMS experiments. For the case of SD interac- 
tions, the contour minimum drops to 5.7xl0 -3 pb, con- 
stituting the most restrictive direct search limit on SD 
WIMP-proton scattering for M w < 60 GeV/c 2 to date, 
and beginning to complement the more sensitive results 
obtained by indirect detection measurements. 

The improved restrictions of the revised Stage 1 con- 
tour are a direct result of the more detailed signal analy- 
sis, improved radio-assays of the shielding materials, and 
the revised nucleation efficiency in the analysis: Stage 
2, with the additional benefit of its improved neutron 
shielding, provides an almost identical sensitivity with 
half the Stage 1 exposure. While the merging may be 
questioned, the results are sufficient motivation for a 
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FIG. 4: various spin-independent contours for Phase II, 
together with those of the leading [l2l Il3t [2l| - [2^ ] spin- 
independent search results; shown are the Stage 2 result, the 
reanalyzed Stage 1 result, and a merging of the two. The par- 
tial contours (a) and (b) are taken from [251 ] and [22], respec- 
tively. The closed areas identified as either CRESST-II [291 ]. 
DAMA/LIBRA [lj or CoGeNT HH represent the regions 
in which possible light mass WIMPS have been respectively 
reported. 



larger exposure measurement with further neutron back- 
ground reduction, and variation of the SDD operating 
temperature/pressure to provide a lower recoil energy 
threshold, towards clarifying the situation. Variation of 
the refrigerant between CaFg, C4F8, CF 3 I or one of the 
other SDDs developed by SIMPLE in recent years @ al- 



lows a variation of detector sensitivities between SI and 
SD sectors [34]. An energy spectrum can in principle be 
obtained, should candidate events be identified, by either 
a temperature or pressure ramping of the SDDs. 

The suggested improvements in this measurement how- 
ever also require a significantly increased active target 
mass in order to be competitive, for which the fabrication 
requirements of the current SDDs become unmanageable: 
a new device, originally prototyped in 2000 |7[ and re- 
prototyped in 2010, based on a large superheated freon 
droplet contained within a gel-sheathed vessel (effectively 
a bubble chamber), is currently completing development; 
it will permit a factor 25 increase in the active mass with 
reduced space requirements. An additional all-around 
60 cm of purified water shielding increases the neutron 
suppression by more than 10 3 , giving the possibility to 
achieve exposures of 10 2 kgd in a few weeks. 
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